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Abstract 

Geometric and electronic properties of folded graphene nanoribbons (EGNRs) 
are investigated by first-principles calculations. These properties are mainly domi¬ 
nated by the competition or cooperation among stacking, curvature and edge effects. 
Eor the zigzag EGNRs, the more stable structures are revealed to be AB stackings, 
while for the armchair types, AA" stackings are more stable. The interlayer inter¬ 
actions and hybridization of four orbitals lead to smaller energy gaps, anti-crossing 
bands, and more band-edge states. Specifically, the broken mirror symmetry in the 
odd-AR stacked zigzag EGNRs is responsible for the spin-up and spin-down splitting 
subbands. All EGNRs are direct-gap semiconductors except that the edge-edge inter¬ 
actions cause the even-AA stacked zigzag EGNRs to exhibit a pair of metallic linear 
bands. The width-dependent energy gaps in the armchair EGNRs can be classified 
into six groups. Eurthermore, there exist rich features in density of states, including 
the form, number, intensity and energy of the special structures. 
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1. Introduction 


In 2004, graphene was first snccessfnlly synthesized (^, and since then has attracted 
considerable attention in the helds of chemistry, materials science and physics [2]-[^. The 
nniqne two-dimensional (2D) hexagonal symmetry indnces rich electronic properties, snch 
as an extremely high mobility [^, an anomalons qnantnm Hall effect [^|^, and ambipolar 
transport phenomena [^. Graphene is viewed as a fntnre electronic material. However, 
its zero-gap property makes actual applications in logic electronics difficult. In order to 
further expand the range of applications, opening an energy gap is critical. Semiconducting 
ID nanoribbons (GNRs) have been the focus of many experimental |9 22 and theoretical 


studies 23 GNRs have been successfully produced from graphene through lithographic 
techniques [9, 10 and from carbon nanotubes by several chemical and physical unzipping 


methods 11-13 . Under different growing environments, GNRs with various structures 


are presented in laboratories, including few-layer 14 , curved 12 , scrolled 15 , and folded 


types 16-^. The unique stacking conhguration and curved structure in folded graphene 
nanoribbons (FGNRs) are worth detailed studies. 

Flat GNRs possess rich geometric and electronic properties, being dominated by their 
edge structures and ribbon widths. Zigzag and armchair GNRs exist according to the achi¬ 
ral boundary condition. A zigzag GNR is a middle-gap semiconductor with two edge-state 
electrons, which have a ferromagnetic arrangement for each edge but an antiferromagnetic 
one across the nanoribbon j^. The zigzag systems are considered as potential materials 
for future applications in spintronic devices, mainly owing to their halfmetallic character¬ 
istics under a transverse electric held [^. On the other hand, all armchair GNRs are 
direct-gap semiconductors, which can be classihed into three types of energy gaps, being 
mainly determined by the number of dimer lines {Ny). The largest (smallest) energy gaps 
are reported to be ai Ny = 3p + 1 {Ny = 3p -|- 2); p is a positive integer. All energy gaps 
are inversely proportional to the ribbon width for all groups [^ . 

In previous research studies, the fundamental properties of GNRs with various cur¬ 


vatures or stackings have not been thoroughly explored. Few-layer GNRs 28 present 
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various stacking configurations. The atomic interactions between two layers strongly de¬ 


pend on the van der Waals interactions. The edge-edge interactions 28,29 are introduced 
by the open edges in some specific configurations. Hence, the electronic properties are 
dominated by both the stacking- and the edge-related interactions. The hybridization of 


four orbitals {2s,2px,2py,2pz) in curved ribbons 24,33 and nanoscrolls causes the complex 
bondings at the inner side of the curved surface. The larger 2pz orbitals on the outer side 
of the curved surface provide a compatible environment for adsorbed atoms, such as H 


and Li 30,31 . As a result, curved systems are promising materials for future hydrogen 
storage or lithium batteries. Due to the cooperation or competition among the stacking, 
curvature, and edge effects, a FGNR with curved and stacking parts (Fig. 1) is expected 
to display intricate and versatile electronic properties. In experiments the synthe¬ 

sized FGNRs have displayed distinct physical characteristics. Experimental measurements 
from scanning tunneling microcopy show that FGNRs can be formed in various stacking 


conhgurations 17,18 . On the other hand, some physical properties of FGNRs have been 


studied in previous theoretical research 32,33 . However, knowledge of the existence of 
more stable FGNR conhgurations, their formation energies, the radius of the curved part, 
band structures, density of states (DOS), and band gaps is still incomplete. Most research 
has focused on the AA stacking, but few investigations have been carried out on the more 
stable AB stackings and other potential conhgurations. 

In this paper, the geometric and electronic properties of zigzag and armchair FGNRs 
are investigated in detail by hrst principles calculations. Many stacking conhgurations are 
considered, not only on the AA and AB stacking conhgurations but also AA! and AA!' 
ones. Among the complex structures, we systematically classify four types of zigzag and 
armchair systems, based on the edge structure, formation energy, ribbon width and stack¬ 
ing conhgurations. The armchair systems are further divided into six groups of width- 
dependent energy gaps on the basis of the stacking conhgurations. Each of these eight 
types of FGNRs presents distinct electronic properties, such as a pair of metallic linear 
bands, metal-semiconductor transitions, splitting of spin-up and spin-down states, and the 
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monotonous width dependence of the energy gaps, etc. Moreover, the predicted main fea¬ 
tures in DOS are validated and compared to experimental scanning tunneling spectroscopy 
(STS) measurements. These rich fundamental features in FGNRs can be expected to pro¬ 
vide potential applications in energy materials, as well as electronic and spintronic devices. 


2. Theory and geometric structure 


We use the Vienna ab initio simulation package 36,37 in the density-functional theory 


(DFT) to investigate the essential properties of FGNRs. The DFT-D2 method 38 is taken 
into account in order to describe the weak van der Waals interactions. The electron-ion 
interactions are obtained by the projector-augmented wave method, and the exchange and 
correlation electron-electron interactions are calculated from the generalized gradient ap¬ 


proximation 39 with Per dew-Burke-Ernzerhof. The wave functions are expanded by the 


plane waves with a maximum kinetic energy of 500 eV. For electronic properties and opti¬ 
mized geometric structures, the hrst Brillouin zone is, respectively, sampled by 300x1x1 
and 12x1x1 k-points in the Monkhorst-Pack scheme. All of the atoms are relaxed until 
the Hellmann-Feynman force is less than 0.03 eV/A. A vacuum space of 15 A is inserted 
between periodic images to simulate an isolated FGNR. 

A FGNR, formed by folding a single flat ribbon, can be viewed as two flat ribbons 
connected by a fractional nanotube in the position to create the curved surface, as shown 
in Figs. 1(a) and 1(f). The open edges are passivated by hydrogen atoms (grey circles), 
and the width is characterized by the number of zigzag or dimer lines. Based on the ribbon 
width, edge structures, and stacking conhguration in the flat region, there are eight types 
of stable FGNRs: (I) AA and (II) AB stackings with even zigzag lines, (III) AA' and (IV) 
AB stackings with odd zigzag lines, (V) AA!' and (VI) AA stackings with even dimer lines, 
and (VII) AA!' and (VIII) AA stackings with odd dimer lines. In the AB stacking, half 
of the atoms share identical (x,y) positions, while all carbons in the AA stacking own the 
same projections. The AA' and AA" stackings are, respectively, obtained by shifting the 
upper layer of the AA stackings along the x-direction with displacements of 6/4 and •\/36/2 
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(a) Zigzag FGNRs (f) Amichair FGNRs 



(c) even-AB (h) even-AA 



(e) odd-AB 


G) odd-AA" 







Figure 1: (a) Side views of even-AA zigzag FGNRs. Top views of zigzag FGNRs with 
respect to (b) even-AA, (c) even-AR, (d) odd-AA', and (e) odd-AR. Side views of (f) 
odd-AA armchair FGNRs. Top views of armchair FGNRs with respect to (g) odd-AA, (h) 
odd-AA", (i) even-AA, and (j) odd-AA". 
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{b is C-C bond length). Furthermore, carbon atoms on both layers have the same chemical 
environment. The hydrogen atoms in the Types (I), (II), (VII) and (VIII) are aligned in 
the same straight line; however, the others types are staggered on different lines. 

3. Optimal geometry and electronic properties 

Each stacking conhguration exhibits unique geometric properties. The curved and flat 
parts of a FGNR, separately, determine the nanotube diameter {Dc) and the interlayer 
distance {Df) (Figs. 1(a) and 1(f)). The former represents the influence of the ribbon cur¬ 
vature, and the latter reflects the stacking effect. Both distances are signihcantly affected by 
the ribbon width and stacking conhguration. As for zigzag FGNRs, the width-dependence 
of Dc and Dj are similar for the four types of stacking conhgurations, as shown in Fig. 
2(a). Dc initially grows with an increase of Ny and then reaches a saturated value at Ny > 
40. For large V^’s, DcS of the even-AA, even-AR, odd-AA' and odd-AR stackings are, 
respectively, 9.2 A, 9.3 A, 8.7 A and 8.5 A. Moreover, D/s are less affected by Ny. This 
weak dependence is rehected in the small variation of average R/’s which are, respectively, 
3.5 A, 3.2 A, 3.3 A and 3.2 A following the same sequence as for the DcS. For armchair 
systems (Fig. 2(b)), on the other hand, the A^^-dependence of Dc is different from the 
above-mentioned relations. Dc weakly depends on the ribbon width except for the even- 
AA stacking, whose Dc decreases with increasing Ny and then reaches a hnal value at Ny > 
40. The DcS of the even-AA, even-AA", odd-AA, and odd-AA" stackings are, respectively, 
7.1 A, 7.0 A, 7.0 A, and 7.2 A. The dependence of Df on the ribbon width is negligible, 
and the R/ of the AA-stacked conhguration is larger than that of the AA"-stacked one, 
regardless of the even or odd width. The perceivable R/’s of the AA and AA" stackings 
are 3.5 A and 3.3 A, respectively. 

When the ribbon width is large enough, both R/ and Dc are nearly constant. R/ is 
dominated by the stacking effect, meaning that even- and odd-AR stackings have shorter 
interlayer distances. This result is in good agreement with those presented in bilayer 
graphene nanoribbons and bilayer graphenes [^. On the other hand, the fact that Dc 
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Figure 2: Width-dependence of and Df for different (a) zigzag FGNRs: even-AA (red 
circles), even-AR (blue triangles), odd-AA' (green squares), odd-AR (orange triangles), 
and (b) armchair FGNRs: even-AA" (red circles), odd-AA (blue triangles), even-AA (green 
squares), odd-AA" (orange triangles). 
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is shorter in armchair FGNRs than in zigzag ones implies that the former are subjected 
to a larger deformation. The deformation, which takes place tangent to the direction of 
the dimer or zigzag line, is responsible for this result. Similar deformations have also been 


observed in zigzag and armchair carbon nanotubes 42 . Another important difference 
related to the curved structure is displayed in the formation energy. 

The formation energy {AE) is dehned as the difference between the total energy of a 
FGNR and that of a flat GNR. This energy is dominated by the interlayer interaction, the 
edge-induced energy and the bending energy, which are mainly contributed by the flat two- 
layer interaction, the edge-edge interaction and the mechanical strain, respectively. The 
third energy is always larger than the sum of the hrst and second ones, so AE is greater 
than zero. Note that the edge-edge interactions are inversely proportional to the edge-edge 
distance (de)- AE decreases with the increasing width, because the edge-edge interaction 
and the bending energy for a saturated Dc are nearly unchanged, while the overlap area 
of two flat sheets grows in a widened ribbon. When the width is sufficiently wide {Ny > 
38 in Fig. 3(a))), AE^s of zigzag systems exhibit the following descending order: even-A A 
stacking > odd-AA' stacking > odd-AR stacking > even-AR stacking. The AR-stacked 
FGNRs have a higher stability since they possess larger interlayer interactions. Moreover, 
the even-AR stacking has a AE lower than the odd-AR one. The d^s of both systems 
are, respectively, 3.1 A and 3.6 A, which leads to the stronger edge-edge interaction in 
the even stacking systems. Such influence is also applied to narrow ribbons; therefore, 
AE associated with the even-AA stackings {de =2.9A) is lower than that of the odd- 
AA' ones {de =3.5 A). On the other hand, the edge-edge interactions do not play an 
important role in the A^y-dependence of AE, but the relationship is dominated by the 
stacking conhgurations; AA"-stacked armchair FGNRs are more stable than AA-stacked 
ones, irrespective of the dimer line number being even or odd (Fig. 3(b)). The edge-edge 
interaction is indistinguishable for the four types of armchair FGNRs with almost the same 
de (~ 3.6A). When the ribbon width is larger, the reason for the armchair FGNRs but not 
the zigzag ones having a higher AE can be explained by the fact that De is shorter in the 
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former than in the latter. 


For zigzag systems, the magnetic properties need to be taken into acconnt. The spin 
arrangements can be changed by the geometric strnctnre. The even-AA and even-AB 
conhgnrations have no magnetic moments. However, the odd-HH' and odd-HH systems. 


like monolayer GNR 23 , present magnetic moments at the open edges. The main reason 
is that the edge-edge interactions can effectively snppress the magnetic moments at a 
snfficiently short distance (<3.5 A 


. In addition, the edge-edge distances for these fonr 


systems are, respectively, 2.9 A, 3.1 A, 3.5 A, and 3.6 A. The magnetic moments cannot 
clearly affect the above-mentioned properties, bnt their effects on electronic properties 
are strong. The systems withont magnetic moments predictably have different electronic 
strnctnres in comparison to monolayer GNR. Moreover, the electronic properties of the 
spin-conhguration-dependent systems can be modnlated by the geometric strnctnre. 

The varions ID FGNRs present the featnre-rich electronic properties. The band strnc¬ 
tnres are signihcantly affected by the edge strnctnre, stacking conhgnration, cnrvatnre 
effect and ribbon width. The zigzag FGNRs exhibit many ID energy bands, as shown in 
Figs. 4(a)-4(d). The occnpied {E'") and nnoccnpied {E^) energy bands are asymmetric to 
the Fermi level {Ep = 0). Most of the snbbands have parabolic energy dispersions. Each 
snbband owns several band-edge states characterized by a local minimnm or maximnm. 
These states are located at = 0 and 1 (in nnits of periodical length), and some are lo¬ 
cated in between them. In the higher-energy region, there are many anti-crossing snbbands 
with band-edge states close to kx =1, mainly owing to the interlayer interactions and the 
enhanced hybridization of the fonr orbitals (2s, 2pa;, 2py, 2pz) on the cnrved snrface. In the 
low-lying snbbands the pair of snbbands nearest to Ep in the regions of 2/3 < kx < 1 and 
0 < kx < 2/3, respectively, possess partially flat and parabolic dispersions. The former 
are mainly contribnted by the local edge atoms (the circle radii representing the contri- 
bntions of the edge atoms). These two snbbands are qnite different among the fonr types 
of zigzag FGNRs. The even-AH stacked systems belong to a ID metal, while the other 
three systems are direct-gap semicondnctors. The valence and condnction bands in the 
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Figure 3: Variations of the width-dependent formation energies for (a) zigzag FGNRs, and 
(b) armchair FGNRs. 
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Figure 4: Band structures of (a) even-AA, (b) even-Ai?, (c) odd-AA'; (d) odd-AB stacked 
zigzag FGNRs. (e) Total charge density of even-AR stacking, (f) Partial charge densities 
with spin arrangements of the odd-AB stacking. The red color indicates spin up and the 


blue color spin down. 
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first systems intersect each other at Ep = 0, and the Fermi-momentnm state is similar 


to the Dirac point in an armchair carbon nanotube 43 . This is because the sufficiently 
strong edge-edge interactions induces more charge carriers which appear in the interlayer 
region near the boundary and thus lead to the loop-like charge density (Fig. 4(e)). For 
the even-AB stackings, the linearly intersecting bands become two anti-crossing parabolic 
bands with a narrow direct energy gap, as shown in Fig. 4(b). The metal-semiconductor 
transition is caused by the variation of the stacking conhguration, which is responsible for 
the fact that a distinct geometric symmetry can open the linearly intersecting bands in the 


same way as in a bilayer graphene nanoribbon 28 and a large closed carbon nanotube 40 


On the other hand, the odd-AA' and odd-AS FGNRs, with magnetism, exhibit two spin- 
dependent band structures (Figs. 4(c) and 4(d)). Both spin-up and spin-down subbands 
are doubly degenerate for the former, while the spin degeneracy is absent in the latter. The 
odd-Ai? systems do not have the mirror symmetry about the z-axis so that the magnetic 
environments are different for edge atoms with spin-up and spin-down conhgurations (Fig. 
4(f)). However, the odd-A A' stackings exhibit the opposite behavior. This can account for 
the spin splitting or degeneracy of the energy subbands. Among the semiconducting zigzag 
FGNRs, the odd-AA' ones have the largest energy gaps. In addition, the energy spacing 
ad kx =1 is smaller in the odd-systems than in the even ones. 

Gompared to the zigzag FGNRs, the armchair ones display an entirely different type 
of band structure. Energy dispersions of the armchair FGNRs are almost independent of 
the stacking conhguration, as shown in Figs. 5(a)-5(h). In the region of < 1 eV, 

each subband exhibits one band-edge state at = 0. However, some of the band-edge 
states are doubly degenerate due to the stacking effect. All the armchair FGNRs belong 
to direct-gap semiconductors, being determined by the = d states. The energy gaps 
are sensitive to changes in the ribbon width and the stacking conhguration. The Ny = 34 
even-AA and Ny = 37 odd-AA stackings possess larger energy gaps. Similar results have 
be predicted for the hat nanoribbon case, where the largest energy gap is from the group 


of Ah = 3p -|- 1; that is attributed to the hnite-width quantum conhnement 23 . When 
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Figure 5: Width-dependent band structures of AA stacked armchair FGNRs: (a)JVy = 34, 
(b)iVj, = 35 , (c) Ny = 36; (d) Ny = 37. The ribbon-dependent band structures of A A" 
stacked armchair FGNRs: (e)W^ = 34, (f)iVj^ = 35, (g)iVy = 36; {h.)Ny = 37. 
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the stacking configuration changes, the Ny = 34 even-AA" and Ny = 37 odd-AA" stacking 
FGNRs also own larger energy gaps. In addition, the AA" stacking has larger gaps than 
those possessed by the AA stacking. 

The versatile width dependence of the energy gaps in various FGNRs deserve closer 
examination. For the zigzag systems. Eg decreases with an increasing ribbon width. The 
energy gaps are strongly dependent on the stacking and curvature effects; therefore, small 
zigzag FGNRs exhibit the following descending order: odd-AA' stacking > odd-AB stacking 
with spin down > odd-AB stacking with spin up > even-AR stacking. The dependence of 
Eg on Ny is weak for the even-AB stackings, while the opposite is true for the odd-AA' 
and odd-AB ones (Fig. 6(a)). This is because the earlier explanation on the formation 
of energy gap in the even-AR stackings considerably differs from what accounts for the 
odd-AA' and odd-AB ones. On the other hand, the energy gaps of armchair FGNRs can 
be further divided into six groups, based on the two types of stacking configurations. R^’s 
corresponding to the AA and AA" stackings, respectively, possess three iVj^-relationships, all 
of which decline with a larger ribbon width as shown in Fig. 6(b). The even-AA stacking 
presents the largest (smallest) energy gap, when the ribbon width is Ny = 6p 4 {6p). 
Moreover, the odd-AA stacking owns a smaller gap in Ny = 6p 5, but the largest energy 
gap lies in Ny = 6p-\-3 or 6p-\-l. Similarly, the even-A A" and odd-AA" stackings are each 
classified into one of three width-dependent energy gap groups. The even-AA" stackings 
belong to the largest energy gap among the discussed armchair system when the ribbon 
width corresponds to the 6p-|-4 group. The smallest energy gap of the even-AA" stackings 
is observed in the 6p -|- 2 group. The odd-AA" stacking for the 6p -|- 1 (6p -|- 3) group owns 
a larger (smaller) energy gap. These complex width-dependences are different from those 
presented in fiat nanoribbons, mainly owing to the combined stacking and curvature effects 
in this unique folding structure. The ribbon width corresponding to the smallest energy 
gap is the 3p -|- 2 group for the fiat nanoribbons; however, this is not exactly the case for 
the FGNRs. Instead, another group associated with the smallest energy gap here is 3p, for 
example, the even-AA and odd-AA" stackings. 
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Figure 6: Variation of width-dependent energy gaps of (a) zigzag, (b) AA stacked armchair 
and (c) AA" stacked armchair FGNRs. 
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The DOS directly reflects the main characteristics of the band structures. The low- 
energy DOS of the even-AA stacked zigzag FGNRs, as shown in Fig. 7(a) by the black 
solid curve, exhibits a plateau near Ep, two prominent peaks (solid and hollow squares), 
and many square-root asymmetric peaks, respectively, coming from the intersecting linear 
bands, the partially rounded bands, and the parabolic bands (Fig. 4(a)). The constant 


plateau indicates the metallic behavior, as seen in a ID armchair carbon nanotube 43 


The two prominent peaks on either side of Ep are mainly contributed by the local edge 
atoms. However, in the even-HH stacking, the continuous plateau of the DOS is destroyed 
and evolves into two separate peaks (red solid and hollow circles in Fig. 7(a)), which 
determine a direct narrow gap. The odd-AA' and odd-HH stackings possess two DOS’s: 
one is associated with spin up (red curve) and the other is related to spin down (blue curve). 
The contributions are not distinguishable to DOS from the spin-up and spin-down states 
for the odd-AA' stackings (Fig. 7(b)), but are distinct for the odd-HH ones (Fig. 7(c)). 
Moreover, regardless of the spin states, the energy spacings of the two adjacent prominent 
peaks (solid and hollow squares) are smaller in the odd-AA' and odd-HH systems than in 
the even-HH ones. The energy gap of the two separate peaks nearest to Ep is the largest 
(smallest) in the odd-AA' stackings {even-AB stackings). The low-lying features in DOS 
are dramatically changed by the edge structure, as displayed for the armchair FGNRs in 
Figs. 7(d)-7(g). The peak energies, numbers and heights near Ep contrast sharply with 
those in the zigzag systems. There is no pair of prominent peaks coming from the edge 
atoms; on the contrary, many peaks originate from the = 0 band-edge states. There 
are some relatively high peaks (green diamonds) are exhibited, mainly owing to the double 
degeneracy. 

The curvature, edge, and stacking effects, which play an important role in the energy, 
number, and height of peaks in DOS, can be verified by the scanning tunneling spectroscopy 
measurements. The tunneling conductance is roughly proportional to the DOS as a mea¬ 
surement reflecting the major structures. The various features of the DOS can be used to 
identify the type of FGNR. The zigzag FGNRs possess a pair of prominent peaks on either 
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Figure 7: Density of states for (a) even-AA and even-^di?, (b) odd-AA', (c) odd-AB stacked 
zigzag FBNRs, and for AA and AA' stacked armchair FGNRs with (d) Ny = 34, (e) 
Ny = 35, (f) Ny = 36; (g) Ny = 37. 
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side of Ep, a characteristic not reflected in the armchair systems. The DOS associated with 
the electron spin is the key to differentiate whether the width of the zigzag systems is even 
or odd. Only the even-AA stacking has free-carrier states at Ep. The splitting of the spin 
degeneracy induces more peaks for the odd-Ai? stackings. For the armchair FGNRs, the 
odd-AA, even-AR", even-RA, and odd-AA" stackings cannot be easily identihed due to the 
similar DOS. However, the energy gaps in DOS are a critical property in conhrming the 
width-dependent six groups. In previous experimental measurements, the curvature effect 


in carbon nanotubes 44 and rippled graphenes 45 and the stacking effect in few-layer 


graphenes 46 have been verihed. As expected, these two combined effects in FGNRs can 


be further examined by the STS measurements. 


4. Conclusion 

The geometric and electronic properties of zigzag and armchair FGNRs are thoroughly 
investigated by calculations based on ab initio functional theory. Eight types of stacking 
conhgurations and various ribbon widths are considered to reveal the fundamental prop¬ 
erties, including optimal geometries, formation energies, charge densities, band structures, 
energy gaps and DOS. FGNRs show many important features, such as the constant Dc 
and Df for larger ribbon widths, the destruction or creation of magnetism, the lower for¬ 
mation energies in the AB or AA" stackings, the metal-semiconductor transitions, and the 
monotonous width-dependence of the energy gap. The predicted results could be verihed 


by experimental measurements of scanning tunneling currents 47 , optical spectra 48 , and 
transport properties (4^. 


The geometric structure of FGNR can be divided into the hat and curved parts; there¬ 
fore, the fundamental properties are strongly dependent on the competition or cooperation 
among the stacking, curvature and edge ehects. The zigzag even-AR stackings or armchair 
AA" stackings are the most stable stacking conhgurations. The electronic structures of 
FGNR in the higher-energy region possess anti-crossing bands because of the interlayer 
interactions and hybridization of four orbitals on the curved surface, while in the lower- 
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energy region the energy gap and energy spacing of the band-edge states are sensitive to 
ribbon width and stacking conhguration. The zigzag even-AA stacking possesses intersect¬ 
ing linear bands and the zigzag even-AB one is a narrow-gap semiconductor. The odd-AA' 
and odd-AB stackings are magnetic materials. The latter has two energy gaps associated 
with the spin-up and spin-down states. All armchair FGNRs are direct-gap semiconduc¬ 
tors. These energy gaps exhibit a monotonic relationship, which decreases with increased 
ribbon width. It is worth noting that the armchair systems can be classified into six groups 
of width-dependent energy gaps. The width-dependencies in the armchair FGNRs are 
different from those in flat nanoribbons. In other words, the smallest energy gap of the 
FGNRs is not necessarily associated with the Ny = ?)p + 2 group. The DOS features, 
including the form, number, intensity and energy of the special structures, directly reflect 
the unique electronic properties. These rich electronic properties and various energy gaps 
are promising for the potential application in electronic devices. 
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FIGURE CAPTIONS 


Fig. 1 - (a) Side views of even-AA zigzag FGNRs. Top views of zigzag FGNRs with 
respect to (b) even-RR, (c) even-RR, (d) odd-RR', and (e) odd-RR. Side views of (f) 
odd-RR armchair FGNRs. Top views of armchair FGNRs with respect to (g) odd-AA, (h) 
odd-AA", (i) even-AA, and (j) odd-AA". 

Fig. 2- Width-dependence of D^. and Df for different (a) zigzag FGNRs: even-AA (red 
circles), even-AR (bine triangles), odd-AA' (green squares), odd-AR (orange triangles), 
and (b) armchair FGNRs: even-AA" (red circles), odd-AA (blue triangles), even-AA (green 
squares), odd-AA" (orange triangles). 

Fig. 3 - Variations of the width-dependent formation energies for (a) zigzag FGNRs, 
and (b) armchair FGNRs. 

Fig. 4 - Band structures of (a) even-AA, (b) even-AR, (c) odd-AA'; (d) odd-AR stacked 
zigzag FGNRs. (e) Total charge density of even-AA stacking, (f) Partial charge densities 
with spin arrangements of the odd-AR stacking. The red color indicates spin up and the 
blue color spin down. 

Fig. 5 - Width-dependent band structures of AA stacked armchair FGNRs: {a,)Ny = 34, 
{h)Ny = 35 , (c) Ny = 36; (d) Ny = 37. The ribbon-dependent band structures of A A" 
stacked armchair FGNRs: {e)Ny = 34, {f)Ny = 35, {g)Ny = 36; (h)iVy = 37. 

Fig. 6 - Variation of width-dependent energy gaps of (a) zigzag, (b) AA stacked armchair 
and (c) AA" stacked armchair FGNRs. 

Fig. 7 - Density of states for (a) even-AA and even-AR, (b) odd-AA', (c) odd-AR 
stacked zigzag FBNRs, and for AA and AA" stacked armchair FGNRs with (d) Ny = 34, 
(e) Ny = 35, (f) Ny = 36; (g) Ny = 37. 
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